Classicality in discrete Wigner functions 
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I. INTRODUCTION 

Continuous- variable Wigner functions W(q,p) have 
been used for a long time to represent quantum systems 
in phase space Q, Q. The Wigner function W(q,p) is 
an alternative complete description of quantum states 
which behaves almost like a phase-space probability den- 
sity. Not only is it real-valued and normalized but it 
also yields the correct value of the probability density 
for the quadrature aQ + bP when integrated along the 
phase-space line aq + bp. However, unlike probability 
densities, the Wigner function can assume negative val- 
ues for some quantum states. This negativity of the 
Wigner function has been considered a defining signature 
of non-classicality (or quantum coherence and interfer- 
ence) 0,0. 

In quantum information science we usually deal with 
systems with a space of states with a finite dimension d. 
For example, for a system of n qubits the dimension of the 
(Hilbert) space of states is d = 2™. For such systems, var- 
ious discrete analogues of the Wigner function have been 
proposed 0, S S II 113 

ITT| and used to investigate 
a variety of interesting problems connected with quan- 
tum computation such as the phase-space representation 
of quantum algorithm s 1121 . separability Qja, quantum 
state tomography [hH Il5lj . teleportation deco- 
herence in quantum walks |Tsj , and error correction 
Here we shall concentrate on a class of discrete Wi gne r 
functions W introduced recently by Gibbons et al. 
This elegant approach seems to be a potentially powerful 
tool to establish connections between phase-space tech- 
niques and problems in quantum information and foun- 
dations of quantum mechanics. 

In this paper we study the set of states with non- 
negative discrete Wigner functions W for all functions 
in the class proposed by jl l| , and the group of unitaries 



that preserve non-negativity of W . Our first result is a 
complete characterization of the set of quantum states 
having non-negative discrete Wigner functions W. This 
is done by proving a conjecture presented by one of us in 
[20| (a related discussion in a somewhat different context, 
usi ng c oncepts in high-dimensional geometry appeared in 
|2lLl22| ). Our proof is elementary and constructive, and 
shows that the only pure states with non-negative W are 
stabilizer states, i.e. simultaneous eigenstates of gener- 
alized Pauli operators [2j|. We then study the group of 
unitaries which preserve non-negativity of W, and prove 
that they form a subgroup of the Clifford group. This 
means such states and unitaries are classical in the sense 
of allowing for an efficient classical simulation scheme us- 
ing the stabilizer formalism. 

The paper is organized as follows. In Section II we re- 
view the discrete Wigner functions W of . In Section 
III we characterize the states with non-negative W, in 
Section IV we discuss positivity-preserving unitary dy- 
namics in phase space and in Section V we summarize 
our results. 



II. DISCRETE WIGNER FUNCTIONS 

In this section we review the class of discrete Wigner 
functions proposed in |ll| and discuss some of their fea- 
tures. 

Let us assume that we are describing a quantum state 
whose Hilbert space dimensionality d is a power of a 
prime number p (d — p n ). In such cases one can introduce 
a phase-space grid with d x d points and label the posi- 
tion and momentum coordinates (q,p) with elements of 
the finite Galois field GF(p n ) 0. At first the use of ele- 
ments of GF{p n ) for both phase-space coordinates could 
be seen as an unnecessary complication, but it turns out 



2 



to be an essential step. The reason is that by doing this 
we can endow the phase-space grid with the same geo- 
metric properties as the ordinary plane. For example, in 
the finite d x d grid we can define lines as solutions to lin- 
ear equations of the form aq + bp = c [where all elements 
and operations in this equation are in GF(p n )). Each 
line will then consist of exactly d points of the grid. The 
field structure of GF(p n ) ensures the validity of proper- 
ties such as: (i) there is only one line joining any given 
two points, (ii) two lines are either parallel (i.e. with no 
points in common) or they intersect at a single point. 
Moreover, it is possible to show that a set of d parallel 
lines (which we will call a sfr , mtzow|ll|') is obtained by 
varying the parameter c in the equation aq + bp — c. Fi- 
nally, the number of different striations turns out to be 
(d + 1). The complete set of (d + 1) striations has been 
studied for a long time in discrete geometry, where it is 
called a finite affine plane p3. |25| . We will label the stri- 
ations with an index k = 1 , . . . , d + 1 and the lines within 
a striation with an index j = 1, . . . , d. In this way the 
j-th line belonging to the K-th striation will be denoted 

as \™. 

A discrete phase space with the above properties was 
used by Gibbons, Hoffman and Wootters in [TlJ to define 
a class of discrete Wigner functions. As mentioned above, 
the crucial property of the continuous Wigner function is 
that its integral along any line A is equal to the expecta- 
tion value of a projection operator P\, i.e. a probability. 
This essential feature is generalized to the discrete case in 
a straightforward way: every line in the d x d phase-space 
grid is associated to a rank one projection operator. As 
noted in this association cannot be arbitrary and 
must obey some simple geometric constraints. For ex- 
ample, we can define a set of d x d unitary operators 
T(q,p) acting on the Hilbert space that faithfully repre- 
sent discrete phase-space translations. For the associa- 
tion between lines and states to respect covariance under 
translations we must impose that the quantum state as- 
sociated to a translated line should be identical to the 
state obtained by acting with the operator T(q,p) on the 
original state. This covariance constraint can be used to 
show the validity of some very significant properties: a) 
the states associated to parallel lines must be orthogonal; 
b) the overlap between states associated to non-parallel 
lines must be equal to l/d. This is important and implies 
that the (d+1) phase-space striations must be associated 
to an equal number of mutually unbiased bases (MUB), 
i.e. bases 

MUBM = {\<t>M),...,\<f>f)} (1) 

such that 

K# ) l# ) )l 2 = ^(i-^ K + ^^- (2) 

As we see, mutually unbiased bases are orthonormal 
bases picked in such a way that any state in one basis is 
an equal-amplitude superposition of all the states of any 
other basis. A complete set of (d + 1) MUB is known to 



exist if the dimensionality of the space of states is a power 
of a prime number. In such case, m any constructions of 
MUB have been proposed jH |$1 M, S El- It has 
been shown that a complete set of (d + 1) MUB for d- 
dimensional systems can be chosen to consist solely of 
stabilizer states, i.e. simultaneous eigenstates of sets of 
(generalized) Pauli operators [23, 123, H3 ■ 

The defining feature of the discrete Wigner functions 
of i s the association between MUB and striations in 
the discrete phase space. As discussed in this can 
be done in a variety of ways and each defines a different 
quantum net [Tl| . which will result in a different defini- 
tion of the discrete Wigner function W. In this paper we 
propose a notion of classicality of quantum states which 
is based on non-negativity of W for all quantum nets ob- 
tainable from a fixed complete set of MUB. It should be 
noted, however, that there have been proposals of crite- 
ria to narrow down the choice of quantum nets: in [l9l | 
the criterion is covariance under the so-called discrete 
squeezing operator; and in the net is chosen so as to 
enforce a natural relation between a separable state's W 
and the W of its subsystems. 

The quantum net is defined by associating each line 

Aj in striation k to a projector Pj = \4>j ){<fij \ onto 
a basis state of basis n. Having fixed a quantum net, the 
discrete Wigner function is uniquely defined by impos- 
ing the condition that the sum of its values along any 
line should be equal to the expectation value of the pro- 
jector corresponding to that line (see ^lj f° r details). 
The resulting Wigner function at any phase-space point 
a = (g, p) can then be shown to be 

W a = Tr (/5i(«)) , (3) 

Ma) = \ £ if'-i), (4) 

where the sum is over projectors associated with all lines 
Aj- containing point a. The construction of the stria- 
tions guarantees that the sum above will contain exactly 
one projector from each basis. 

The operators A(a) are known as phase-space point 
operators and form a complete basis for the space of op- 
erators, which is orthogonal in the Schmidt inner prod- 
uct (i.e. Tr (^A(a)A(fi)j = S a ^/d). We can rewrite the 

expression for the Wigner function at phase-space point 
a using the probabilities associated with the projectors 
p(«0. 

3 ' 

pf=Tr(pPf). (5) 

In terms of these probabilities, the Wigner function at 
the point a takes the form 

W * = \\ £ Pi " 1 ]- ^ 
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The discrete Wigner function W can be shown to have 
many of the features of the continuous Wigner function 
W(q,p) [lU: it is real (but can be negative), normalized, 
and its values are obtained through eq. © from measure- 
ments onto MUB. Here the MUB projectors play the role 
that the quadratures aQ + bP play in W(q,p), forming 
a particularly symmetric set of observables whose mea- 
surement results completely characterize the state (in a 
process known as quantum tomograph y). F or a discussion 
of further properties of W see [ll|, [ljj, |2i| • 

In the discussion that follows we will often be repre- 
senting quantum states using the probabilities pj . As 

the projectors P.- form an over-complete basis for the 
space of density matrices, these probabilities completely 
characterize the state. Since for any striation y\ p^ — 
1, there are only (d— 1) independent probabilities for each 
basis, resulting in a total of (d— l)-(d+l) = (d 2 — 1) inde- 
pendent probabilities, exactly the number of real param- 
eters necessary to describe a general normalized mixed 
quantum state in <i-dimensional Hilbert space. Each 
quantum state is represented by a point p in this (d 2 — 1)- 
dimcnsional probability space. 

As mentioned above, for power-of-prime d it is possible 
to build a complete set of (d + 1) MUB using only sta- 
bilizer states, i.e. joint eigenstates of generalized Pauli 
operators. Let us discuss more explicitly such construc- 
tions for the case d = 2 n , i.e. n qubits (see for more 
details). In order to define a complete set of (2" + l) MUB 
we start by partitioning the (4 n — 1) Pauli operators (ex- 
cluding the identity) into (2™ + l) sets S { of (2™-l) Pauli 
operators each. We will require that the Pauli operators 
in each set Si be mutually commuting, but otherwise the 
partitioning can be completely arbitrary. If we add the 
identity and a ±1 phase to the Pauli operators in each 
set Si, each will form a maximal Abelian subgroup of the 
Pauli group. The joint eigenstates of each such set 5, 
form a basis for the Hilbert space, and due to properties 
of the Pauli operators the (2™ + 1) bases thus defined can 
be shown to be mutually unbiased [28| . 

The phase-space construction provides a natural proce- 
dure for partitioning the Pauli group into disjoint, mutu- 
ally commuting sets. The idea, which is worth reviewing 
here, was described in ^2 and further elaborated in [lflj . 
Pauli operators represent phase-space translations and 
can be labelled using binary n-tuples p and q (n-tuplcs 
q and p contain the coordinates of the field elements q 
and p in a given basis as described below). Each Pauli 
operator can be written as 

n-1 

f{q,p) = J{X?Z?e i ^, (7) 

i=0 

where Xi and Zi stand for the Pauli operators on qubit i, 
and the phase is chosen so as to make the operators Her- 
mitian. The definition above will be written in shorthand 
as 

f(q,p)=X i 'ZPe^^P. (8) 



The condition for two Pauli operators to commute turns 
out to be 

[f(q,p),f(q',p')]=OiSq-p'-p-q' = (mod 2). (9) 
Let us consider a set of (d — 1) Pauli operators 

S (S,b) = {f(aM j ,bM j ) 3 = 0, 1, • • • , d - 2} , (10) 

where M is an arbitrary binary matrix, M is its transpose 
and a, b are binary n-tuples. Any two operators of this 
set commute. It is interesting to note that S^ s a forms a 
maximal Abelian subgroup of the Pauli group if and only 
if M is a generating element of the matrix representation 
of the field GF(2 n ). This can be seen as follows: the 
product of two elements of S,- a is itself an element of 

this set (up to a sign) iff the matrix M is such that for 
every pair of integers j,f, there is a third integer j" 
such that AP + M j ' = M j " [where < < {d - 

2) and j 7^ j']. Moreover, for the set to have exactly 
(d — 1) different elements, the matrix M should be such 
that all powers M 3 for j = 0, . . . , (d — 2) are nonzero 
and different from each other. For M satisfying these 
conditions, it can be seen that M d_1 = 1. Therefore, the 
elements of the set {0, 1, M, M 2 , . . . , M d ~ 2 } form a finite 
field, and we see that the matrix M and its powers form 
a matrix representation of GF(2 n ). A possible choice for 
M, used in [l^, is the so-called "companion matrix" of 
the primitive polynomial which defines the product rule 
in the field. With such a matrix we can build (d + 1) 
disjoint sets of commuting Paulis of the form S, s a by 

choosing the binary n-tuples (a, b) as explained below. 

The association between each phase-space point (q,p) 
and a Pauli operator T(q,p) must respect the covariancc 
of the construction under phase space translations. This 
is done as follows: the line formed by all phase-space 
points satisfying the equation bq + ap = c is invariant 
under phase-space translations of the form q' = q + ato 3 , 
p' = p+buj 1 (where lu is a generating element of the field). 
To this phase-space translation we must associate an op- 
erator acting in Hilbert space. The natural identification 
is to associate this with the operator T(aM 3 ,bM 3 ). Here, 
the choice of n-tuples a and b is arbitrary. The important 
point is that once this choice is made [i.e., once we arbi- 
trarily assign two n-tuples to the point (a, b)] we repeat- 
edly apply the matrix M (M) to the position (momen- 
tum) coordinates to obtain the n-tuples parametrizing 
the Pauli operators associated to the other phase-space 
points [Tflj . In summary, this construction associates an 
operator T(q,p) to every phase space point (q,p) in such 
a way that the elements of the Abelian subgroup S,~ a 
are associated to points in phase space that belong to 
the ray defined by the equation bq + ap = (a ray is 
defined as a line that contains the origin). In 19] it 
was shown that by varying the n-tuples (a, b) one can 
construct only (d + 1) different sets S (s ^. If we define 
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the two n-tuples 1 = (1000 ■ ■ • 0) and = (0000 • • ■ 0), 
these maximal mutually commuting sets of Pauli opera- 
tors can be conveniently built by choosing (a, b) = (1,0) 
(which we will associate with the horizontal striation), 

(a,b) = (0,1) (the vertical striation) and (a,b) = (1,6) 
for b ^ (the other striations). 

Thus, the mapping between striations and MUB is 
naturally determined by the phase-space construction, 
as lines which are invariant under the transformations 
q' = q + auji , p 1 = p + buji must be associated to states 
which are invariant under the corresponding transforma- 
tions in Hilbert space, that is, the translation operators in 
S^ s g-j . Therefore, the lines of the form bq+ap = c must be 
associated to common eigenstates of the set S^ s ^ . How- 
ever, there is no criterion telling us how to associate each 
line in a striation with a projector in the corresponding 
basis. We can count the number of possible quantum 
nets as follows: for the ray of a given striation there are 
d possible projectors to choose from; once this choice has 
been made the condition of covariance under translations 
determines which projector should be associated to each 
of the other lines in the same striation. As there are 
(d+ 1) rays, the number of possible associations between 
lines and projectors is d d+1 , each of which defines a dif- 
ferent quantum net, leading to a different definition of 
the Wigner function. The projectors associated to each 
of the lines in the vertical and horizontal striations can 
be chosen in such a way that the coordinates of each line 
correspond to the eigenvalues of the Paulis generating 
the set (the single qubit Paulis Z and X, respectively). 
Then, there are still possible quantum nets, each 

of them given by a particular choice of projectors to be 
associated to the rays of the remaining oblique striations. 

There is a closely related methodology for constructing 
the Wigner functions for general dimension d = p n which 
emphasizes the link between the exponents of the Z and 
X operators and the finite geometry of V2[GF(j> n )], the 
two-dimensional vector space over the field GF(p n ). One 
defines lines, rays and striations in this two dimensional 
space and then, using the properties of the algebraic field 
extension, defines an isomorphism with V2 n [GF(p)]. Vec- 
tors in this second space serve as exponents of the Z and 
X operators, and the commuting classes of generalized 
Pauli matrices correspond precisely to parallel lines in 
a striation in the first vector space. Details of this ap- 
proach and a methodology for assigning projections to 
lines are given in |30| . 

III. STATES WITH NON-NEGATIVE WIGNER 
FUNCTIONS W 

Following let us now characterize the set of states 
having non-negative discrete Wigner functions W simul- 
taneously in all definitions proposed by Gibbons et al. 
|ll| for power-of-prime dimension d. 

Definition: The set Cd is defined as the set of (pure 



or mixed) density matrices of systems in a d- 
dimcnsional Hilbert space having non-negative dis- 
crete Wigner function W in all phase-space points 
and for all definitions of W using a fixed set of 
mutually unbiased bases. 

By definition, the set Cd is specified as the intersection 
of a number of half-spaces in the (d 2 — l)-dimensional 
p-space. From JBJ it can be seen that each half-space 
inequality is of the form 

E ^ * i. en) 

where the probabilities appearing in the sum are associ- 
ated with the lines containing phase-space point a, and 
hence depend on a and on the quantum net chosen. The 
intersection of the half-spaces defined by these inequal- 
ities is a convex polytope in p-space, given in an H- 
description (H standing for "Half-space"). Any convex 
polytope also admits an alternative V-description (V for 
"Vertices" ) , consisting of the list of vertices whose convex 
hull defines the polytope. 

Galvao showed that for d < 5, the H-polytope Cd has 
a V-description whose vertices are the MUB projectors 
poj . and conjectured this would also be true for general 
power-of-prime d. A geometrical argument showing the 
validity of this conjecture was given in pH |2^. Let us 
now provide a constructive, analytical proof. 

Theorem 1 For any power-of-prime Hilbert space di- 
mension d, the H-polytope Cd is equivalent to the V- 
polytope C v having the MUB projectors as vertices. 

Proof: 

Let us prove the theorem by first showing that the 
V-polytope C v is contained in Cd, and then the con- 
verse. From j 1 1 1 ] we know that the Wigner function for 

any MUB projector pj K ' = \<t>j K ^}{<j>j K ^\ is non- negative. 
Since the Wigner function depends linearly on the den- 
sity matrices, W is non-negative also for any state in the 
convex hull of the Pj^ . This shows that any state in C v 
is also in Cd, as we wanted to prove. 

Let us now prove the converse, i.e. that polytope C v 
is contained in polytope Cd- What is required now is to 
show that any state p 6 Cd can be written as a convex 
combination of the projectors Pj ■ 

p = EE4 K) ^ d 2 ) 

K=l j = l 

(k) 

with all Cj > 0. Note that the decomposition is not 
unique. 

Let us start by considering the general expression for 
Wigner function W at phase-space point a [eq. ©]. 
Given a state p there is some Wigner function defini- 
tion which, at some point a, evaluates to a minimum 
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value among all definitions and all points a. This hap- 
pens when the expression for W a is such that the sum 

(|fi*|) includes only the smallest probability pj from each 
MUB k. Let us denote these VF-minimizing probabili- 
ties p* = minj-jj^}. States in Cd are those for which 
all expressions of the form 10 are non-negative, and this 
happens if and only if the expression for W involving 
only the pi K ^ is non-negative. In other words, a state has 
non-negative W in all definitions if and only if: 

d+l 



Eq. H18|) tells us that each coefficient can be writ- 
ten as the sum of a non-negative term (p^ — P*^^ plus 

a (possibly negative) constant x^ K K We can show, how- 
ever, that the sum of those constants x^ K > has to be non- 
negative. We do that by using the normalization condi- 
tion on eq. JTSJl: 



£ 



» 



E^-E^+E 



i 



« -3 



£^ K) > i. 



(13) 



£=1 



This is our hypothesis. 

Any density matrix p can be expanded in terms of the 



projectors P, 



as in eq. (|12|) . with real 



(but possibly negative) coefficients cj- . A first constraint 

on the coefficients ci K ^ comes from the requirement that 
Tr(j5) = 1. Using property J2J of the MUB we can com- 
pute the trace, obtaining 



d+i-dj2p ( * K) +dJ2 x(K) = 1 

K K 

=► EX 10 = 5> w - 1. (2°) 



Now remember that our hypothesis is that Y^ K P* — 1> 
which implies that ^ K > 0. Let us now use this fact 
and expression l|18|l to obtain an expansion of the den- 
sity matrix p in terms of the projection operators Pj K \ 
but now with non-negative coefficients only. Plugging eq. 
(|T%|) into eq. (jT2|l we obtain: 



Tr(p) 



d+l d 

EE 

K=l j=l 



cf = 1. 



(14) 



(21) 



Now let us use eq. © to calculate explicitly from Usin S thc fact that Ej Pj = we can rewrite this as 
eq. (|12|) , so as to obtain relations between the coefficients 



c^f and the probabilities p- • 



» 



3 ' 

d+l d 



E 



^WhEE* 

/i— 1 m — 1 
i2 



(M) 



E E (#° > I" + E ^ I ^ i« 



p 



(«) 



(22) 



EE 



c O)± 4. r 
m d 3 



(«) 



(15) 



where we have used the condition of mutual unbiasedness 
of the bases [eq. (J2J)]. Now we can use the trace condition 
(114(1 to rewrite this as 



P 



O) 



or 



» 



„(«) 



» 



1 

d 

1 

d 



(16) 



(17) 



Let us add = p* — p*" -1 to the right-hand side of the 
equation above, to obtain 



» 



» 



Ik) (re) 
• - Pi' 



,0) 



with 



-Veto 
d^ m ' 



(18) 



(19) 



where we defined x = ^2 K x < - K \ Note that p^ > pi^ 
by definition, and our hypothesis guarantees that x > 0. 
What we have now is then an expansion of p in terms of 
the MUB projectors using only non-negative coefficients. 
QED 

Our proof above is constructive - for any state in Cd 
we can use equations (jHJ and 1)20(1 to obtain a convex de- 
composition of the state in terms of the MUB projectors 
and their associated probabilities, given by eq. ((22(1 . As 
noted in some non-physical states (i.e. described by 
non-positive Hermitian matrices) can have non-negative 
Wigner functions in a single definition of W. Theorem 
n shows that imposing non-negativity of W for all def- 
initions of W is sufficient to guarantee that the set Cd 
contains only physical states. 

In the light of Theorem ^ above, let us now discuss in 
which senses states with non-negative Wigner functions 
are classical. We have defined the set Cd of states of 
a d-dimensional system with non-negative Wigner func- 
tion W in all definitions. Theorem^Jproves that the only 
pure states in Cd are the MUB projectors, which can al- 
ways be chosen to be stabilizer states, i.e. simultaneous 
eigenstates of Pauli operators [23, |2!J |3(| • The stabilizer 
formalism then provides us with a way to represent pure 
states in Cd using a number of bits which is polynomial 
in the number of qubits [3l| . This contrasts with gen- 
eral quantum states whose classical description requires 
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an exponential number of bits. We thus see that states 
whose discrete phase-space description avoids Feynman's 
"negative probabilities" Q are classical also in the sense 
of having a classical-like short description. 

Classicality witnesses are observables whose expecta- 
tion values, if negative, indicate some non-classical prop- 
erty such as entanglement (see [j^ H^)- Our Theorem 
gives such an interpretation to the phase-space point op- 
erators A(a): negative expectation values indicate (i.e. 
witness) non-classicality in the sense we discussed above. 

The stabilizer formalism provides us with a framework 
in terms of which pure states in the set Cd have an effi- 
cient classical description. Other choices of frameworks 
are possible, each choice resulting in a different set of 
quantum states with efficient classical descriptions. One 
example are the (mixed) separable density matrices of a 
collection of qubits, each of which has an efficient classi- 
cal description in terms of single-qubit pure states. An 
efficient description, however, does not guarantee the ex- 
istence of an efficient simulation scheme for the dynamics; 
the dynamics of separable mixed states in NMR quantum 
computation experiments provides us with an example of 
this problem PJ| . In the next section we discuss the issue 
of simulability of unitary dynamics within our set Cd- 

Given these observations, it is not surprising that pure 
states in Cd can behave non-classically in other ways, 
that is, with respect to other frameworks. For example, 
states in Cd can be highly entangled, allowing for proofs 
of quantum non-locality and contextuality. 



IV. UNITARIES PRESERVING 
NON-NEGATIVITY OF W 

In continuous phase-space we can define classical uni- 
taries as the group of unitaries which preserve non- 
negativity of the Wigner function W(q,p). It has been 
shown that this group is formed by all unitaries gener- 
ated byHamiltonians which are quadratic forms in phase 
space 0| • In this section we obtain an analogous result 
for our discrete Wigner functions W: using the 'classical' 
pure states in Cd, we define and characterize the group of 
unitaries {U c } that map pure states in Cd to other pure 
states in Cd- In other words, we characterize the group of 
'classical' unitaries {U c } that preserve non-negativity of 
W for all quantum nets obtained from a fixed complete 
set of MUB. 

The structure of the group {U c } may depend on the 
particular complete set of MUB we choose to define our 
discrete Wigner functions. We prove that for any MUB 
construction using Pauli operators, the group {U c } is a 
subgroup of the Clifford group (the group of unitaries 
mapping Pauli operators to Pauli operators under conju- 
gation pi)). For the particular construction in [lll[l9) we 
present some unitaries in {U c } and discuss their action 
in phase space. 



A. {U c } is a subgroup of the Clifford group 

Let us consider the (4" — 1) Pauli operators acting 
on the 2™-dimensional Hilbert space of n qubits (exclud- 
ing the identity). We can partition these (4™ — 1) Pauli 
operators into (2" + 1) sets Si of (2™ — 1) commuting 
Paulis each. The joint eigenstates of the sets Si form a 
complete set of (2™ + 1) MUB, as discussed in section 
ITT1 We can partition the Paulis in many different ways; 
each such partitioning defines a different complete set of 
MUB, which will be denoted as Bj. In this section we 
show that the 'classical' unitaries {U c } mapping MUB 
in a partition Bi to MUB in the same partition form a 
subgroup of the Clifford group. 

The strategy is as follows. We will consider a slightly 
more general problem, which is to characterize unitaries 
mapping MUB defined by an arbitrary partition B\ of 
Pauli operators to MUB defined by a second partition 
£?2- A general unitary U mapping B\ to B2 will in 
particular map two bases in B\ to two other bases in 

B2- Let us name them (Si G Bx) & (S 2 G B 2 ), 

(Ti 6 Bi) ^> (T 2 6 5 2 ). The first step is to prove there 
are two Clifford unitaries Cj(j = 1,2) that map basis Sj 
to the computational (Z) basis, while mapping basis Tj 
to the X-basis. These standard Clifford unitaries are the 
key to the proof. This is because U is Clifford if and only 
if U = C2UCI is Clifford. So it is enough to show U is 
Clifford (done in Theorem^, which is easier as by con- 
struction U are unitaries that preserve both the Z basis 
and the X basis. 

With this more general result in hand, we can consider 
the case when the two partitions are one and the same 
(Bi = B2), and we will have what we wanted to prove, 
i.e. that our 'classical' unitaries {U c } are Clifford group 
operators. 

We wish to show 

Theorem 2 Let U be a unitary transformation that 
maps a complete set of Pauli MUB B\ to a second com- 
plete set of Pauli MUB B 2 . Then U is in the Clifford 
group, up to a global phase. 

The first step involves proving the following Lemma: 

Lemma 1 Let S and T be two maximal Abelian sub- 
groups of the Pauli group, with S PI T = {1}. Then there 
exists a Clifford operation which maps S 1— > {Z\, . . . , Z n ), 
T h-> (Xl, . . . ,X n ). 

Proof: 

Since S and T are maximal Abelian subgroups with 
trivial intersection, it follows that no (non-identity) ele- 
ment of T commutes with every element of S (or vice- 
versa). Let {Mi i = l,...,fi} be a set of generators 
of S. For any particular clement N G T, we can define 
the syndrome <j(N) which is an n-tuple whose i-th com- 
ponent is given by Oi = c(N,Mi) i = 1, ...,n. Here 
c(N,M) = if N and M commute and c(N,M) = 1 
if N and M anticommute. Then it follows that if 
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N,N' £ T, N ^ N', then a(N) ^ a(N') (since oth- 
erwise a(NN') = a(N) + a(N') = 0, and NN' would 
commute with every element of S). 

In particular, since there are 2™ elements of T and 2™ 
different possible values of a, it follows that each value 
of a is used exactly once. Thus, we can choose Ni G T 
such that a(Ni) — el (where el is the vector that is 1 
in the i-th position and elsewhere). That is, Ni anti- 
commutes with Mi and commutes with Mj (i ^ j). The 
A/j's are independent (because their a vectors are inde- 
pendent) and they commute with each other (because T 
is Abelian). Therefore, the set of Mi's and iVj's have the 
same commutation/anticommutation relationships as the 
Zi's and the X^s, so there exists a Clifford group oper- 
ation that maps Mi i— > Zi and Ni i— > Xi. This provides 
the appropriate map on S and T. QED 

This Lemma can be adapted so it applies also to d- 
dimensional registers, the main difference being that the 
syndrome function a(N) assumes values which are vec- 
tors modulo d (see |35|)- 

An immediate consequence of this lemma is that for 
any complete set of Pauli MUB B, we can choose any two 
of its bases, represented by stabilizers S and T, and find a 
Clifford group operation that will map B to another com- 
plete set of MUB containing the bases (Zi, . . . , Z n ) and 
(Xi, . . . , X n ); and in particular, this Clifford group oper- 
ation will map S to (Z\, . . . , Z n ) and T to (Xi, . . . , X n ). 

Therefore, if we have a general unitary U that maps 
Pauli MUB B\ to Pauli MUB B 2 , we can choose bases 
St,Ti G Bt with S 2 = U(Si), T 2 = U(Ti) (so S 2 ,T 2 G 
B 2 ), and then find Clifford operations C\ and C 2 which 
map C\ : Si h-> (Zi,...,Z n ), C\ : T x i-> (Xi,...,X n ), 
C 2 : S 2 ^ {Zx, . . . , Z n ), C 2 :T 2 ^ (X x , ...,X n ). Then 
it follows that C 2 UC\ : (Z 1 ,...,Z n ) ^ (Z u ...,Z n ), 
C 2 UC\ : {X x , . . . , X n ) ^ (X u . . . , X n ). Denoting U = 
C 2 UC\, it is easy to see that U is a Clifford group oper- 
ator iff U is a Clifford group operation. Thus, to prove 
Theorem 12 it will be sufficient to prove 

Theorem 3 If U is a unitary operation which preserves 
both the Z basis and the X basis (i.e., maps eigenstates of 
(Z\, .... Z n ) to other eigenstates of this set of operators, 
and the same is valid for eigenstates of (X\, . . . ,X n )), 
then U is a Clifford group operation, up to a global phase. 

Proof: 

Since U preserves the Z basis, it has the form of a 
classical gate with possibly some phases changed: 



Therefore, 



(5) 



(23) 



with g(z) a permutation of the 2™ possible values of z. 

In terms of the Z basis, we can expand elements of the 
X basis as follows: 



\x) 



X 



t(x-z) 



(24) 



U\x) x = J2< 



= y^ e i(TX-g (z)+4>{g ( 2 5) 
z 

In order to preserve the X basis, we need 

U\S) X = e l6 ^\h(x)) x = ^e l( - e( - s)+ ^ {s) - S) \z) z , (26) 

z 

where h(x) is a permutation of the values of x. Equating 
(ESI and (EH), we find 

■nx-g _1 (z) + 4>(g ~ 1 (z)) = 6{x) + irh(x) ■ z. (27) 

This must be true for all x and z. Plugging in x — 0, we 
find 

cf>(g- 1 (z))=9 + 7rh a -z, (28) 
where 8q = 8(0) and ho = h(0). Therefore 



ttx ■ g 1 (z) = 7r h(x) — ha ■ z + [0(x) — 9q 



(29) 



Of course, eqs. are understood to be modulo 2ir. 

It then follows that g ~ 1 {z) must be affine in z: 



g 1 (z) = Az + b, 



(30) 



(where A is an invertible n x n binary matrix) and thus, 
by H28|) . 4>(z) is also affine in z: 



4>(z) = 7TC • Z + d, 

with Ac = ho and irb ■ c + d = 9q. 
Thus we find 



U\z) z = e l(7TC - 



+d ^\A- 1 z-A- 1 b). 



(31) 



(32) 



We can easily identify this as a Clifford group opera- 
tion, up to the global phase e ld : \z) 1— > can be 
performed with CNOT gates, \z) t—> \z — A~ 1 b) can be 
performed with X operations, and \z) >— > e l7r ( c ' z ) |^) can 
be performed with Z operations. QED 

The proof for d-dimensional registers is almost identi- 
cal, except that we must replace 7r everywhere with 2n/d, 
and we need scalar multiplication gates as well as SUM 
gates to perform \z) >— ► |^4 _1 5) [3^. 

For prime Hilbert space dimensions there is a unique 
Pauli MUB construction that uses all stabilizer states 
p9). In that case our set of pure classical states coin- 
cides with the set of stabilizer states, and our group of 
classical unitaries {U c } coincides with the Clifford group. 
This is not the case for power-of-prime dimensions, where 
a complete set of Pauli MUB contains only a proper sub- 
set of the stabilizer states, resulting in classical unitaries 
{U c } which form a proper subgroup of the Clifford group. 
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The Gottesman-Knill theorem [23( states that Clifford 
group operations on stabilizer states can be simulated ef- 
ficiently on a classical computer. Thus, our Theorems 
n and [21 guarantee that the group of classical unitaries 
{U c } applied on the set of classical pure states in C d can 
be efficiently simulated, i.e. with a number of time-steps 
that increases only polynomially in the number of qubits. 
This is to be contrasted with general quantum computa- 
tion, which uses states and operations outside of our clas- 
sical sets, and which is thought to provide exponential 
speedup for some problems. The necessity of negativity 
of W for achieving universal quantum computation had 
been noted in pp| for a particular computational model 
proposed recently by Bravyi and Kitaev • 

Our results point to an interesting convergence be- 
tween two different notions of classicality. The first de- 
fines classical states as those whose description can be 
made in terms of non-negative quasi-probability distribu- 
tions, in this case the discrete Wigner functions of [ll|. 
The second is motivated by quantum computation: clas- 
sical states and operations are those which can be effi- 
ciently simulated on a classical computer. For a related 
discussion of simulability in the context of continuous 
variables see |37| . 

In this section and the previous one we only made 
claims of classicality for pure states in Cd and their asso- 
ciated unitary dynamics. The problem seems to become 
much more involved when we consider mixed states in 
Cd and their associated dynamics, which in this case will 
be (in general non-unitary) completely positive maps. It 
is not clear whether an efficient simulation scheme for 
this more general definition of classical dynamics can be 
devised. 



B. Unitaries in {U c } and their action in phase space 

We have just shown that when we build a complete set 
of MUB using Pauli operators, the 'classical' unitaries in 
{U c } turn out to form a subgroup of the Clifford group. 
The exact characterization of this subgroup will depend 
on which Pauli MUB construction we pick. In this sec- 
tion we restrict ourselves to the construction for N qubits 
sketched in section |HJ and present some 'classical' uni- 
taries in {U c } together with their associated action in 
phase space. 



1. Discrete phase-space translation operators 

In ^3 it was shown that the Pauli operators act as dis- 
crete phase-space translations mapping phase-space lines 
into other lines. This means that the Pauli operators 
themselves are in our group {U c }. Translation operators 
T(g,p) operate on quantum states in such a way that 
their Wigner functions are transformed as flows: each 
phase-space point operator is mapped into another one 







V 





























FIG. 1: Operators in {U c } map translations T(q,p) into 
T(q',p') (up to a sign); representation of the transformation 
of translations under the action of: a) translation operators; 
b) discrete squeezing; c) Fourier transform. 



because of the covariance condition imposed on the quan- 
tum net. Thus, the effect of a translation operator on a 
state's Wigner function is a translation in phase space, 
and in this sense its action is "classical-like" . Since trans- 
lation operators act as flows in phase space, they preserve 
positivity of W for any single association between lines 
and MUB projectors. 



2. Discrete squeezing operator 

The discrete squeezing operator 19] maps the horizon- 
tal and the vertical striations into themselves, while cy- 
cling through all oblique striations. An explicit Clifford 
circuit for U s is given in 0] . When acting on translation 
operators, U s maps them into other translations in a way 
that resembles a squeezing flow in phase space (see Fig. 
l.b): 



U s f{q,p)Ut = ±T{qM,pM- 1 ). 



(33) 



Besides covariance with respect to the discrete phase- 
space translations, we can impose on the quantum net 
also the constraint of covariance under U s . In doing so, 
the freedom in picking the quantum net will be limited to 
the choice of which MUB projector to associate to a fixed 
oblique line, since the covariance requirement determines 
all other associations. In this way, the number of possible 
choices is greatly reduced from d d_1 to d. 

If we choose a quantum net which is covariant under 
the squeezing operator, it can be shown that U s will map 
phase-space point operators into other point operators, 
i.e. U s acts like a phase-space flow. This may not be 
the case if the quantum net is not chosen to be covariant 
under U s . 

By definition, the group {U c } consists of unitaries that 
preserve non-negativity of the Wigner function W for all 
possible quantum nets. This does not imply that opera- 
tors in {U c } will preserve positivity for any single defini- 
tion of W. This is because some states may have positive 
W for a single definition, but negative W for other def- 
initions (and hence lie outside the set Cd)- In the case 
of U s , preservation of positivity for each single definition 
of W is only guaranteed when U s acts as a flow in phase 
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space, and this only happens when the quantum net is 
chosen to be covariant with respect to U s . 

3. Finite Fourier transform 

The finite Fourier transform F [3^ maps the horizon- 
tal and the vertical striations into one another; oblique 
striations are interchanged in pairs, and one of them [the 
"main diagonal" , which corresponds to the eigenstates 
of the set of translations obtained by setting a = b in 
eq. ljTU)l] is mapped into itself. For the particular case in 
which the canonical basis of the Galois field is self-dual, 
F is just the Hadamard transform. The effect of F on 
the translation operators is - up to a sign - a reflection 
with respect to the main diagonal of phase space (see Fig. 
l.c): 

FT (lM j , lM fc )Ft = ±f(lM k , lM j ) (34) 

The fact that F interchanges translation operators by 
a reflection might suggest that its action on the states 
could be analogous, that is, that F could reflect a state's 
Wigner function with respect to the main diagonal, per- 
haps for some particular quantum nets (as is the case 
with U s ). 

For F to act on lines as a reflection with respect to 
the main diagonal, there should be one MUB projector 
associated to the axis of reflection, and F should map 
this projector into itself. This projector should, then, 
be a common eigenstate of F and all the Paulis that 
define the basis to which the state belongs. Using the fact 
that F anti-commutes with some of them, and that the 
eigenvalues of F and the Paulis are different from zero, 
it can be seen that no state can fulfill this requirement. 
Thus, there is no association between lines and MUB 
projectors that makes the action of F on the Wigner 
function be a reflection flow. 

Moreover, it can be seen that there is no quantum net 
for which F acts as a flow in any way, because for F to 
be a flow it should map phase-space point operators into 
other point operators. For this to happen, the (d + 1) 
lines that intersect in any given point must be mapped 
by F into other lines that intersect in only one point. The 
vertical and the horizontal rays (i.e. lines containing the 
origin) are interchanged by F, so the other rays must 
be mapped into rays too (so that all the resulting lines 
intersect at the origin) . This requires the ray in the main 
diagonal to be mapped into itself, and, as pointed out in 
the previous paragraph, this cannot be achieved. 



Therefore, F provides an example of an operator in 
{U c } which cannot be interpreted in terms of a flow for 
any choice of associations between lines and states, and 
so has no obvious continuous phase-space analogue. 

V. CONCLUSION 

We have characterized the set Cd of states whose dis- 
crete Wigner functions W (as defined in [llj) are non- 
negative. We showed that the only pure states in Cd 
are the mutually unbiased bases projectors used to de- 
fine W, as conjectured in [2(|. Since these projectors 
can always be chosen to be stabilizer states, they admit 
an efficient classical description using the stabilizer for- 
malism. Moreover, we proved that the unitaries which 
preserve non-negativity of W for all such functions W 
form a subgroup of the Clifford group. It is known that 
Clifford operations on stabilizer states can be simulated 
efficiently on a classical computer. We have thus iden- 
tified a relation between two different notions of classi- 
cality: states which are classical in the sense of having 
non-negative quasi-probability distributions (the discrete 
Wigner functions of (T^) can also be simulated efficiently 
on classical computers. Since general quantum compu- 
tation is thought to be hard to simulate classically, our 
results mean that negativity of W is necessary for expo- 
nential computational speedup with pure states. 

There are many open problems worth investigating. 
The complete characterization of non-negativity preserv- 
ing unitaries for different constructions of complete sets 
of MUB is still unsolved. It would also be interesting if 
one could relate non-classicality to negativity of If in a 
quantitative way. Another research direction is to inves- 
tigate the relationship between W and a notable open 
problem, that of the existence of complete sets of mutu- 
ally unbiased bases for general Hilbert space dimensions 
(see 0,H3,|2^|). The original idea behind continuous- 
variable Wigner functions was to help visualize quantum 
dynamics in the familiar framework of classical phase 
space. Some research has been done on the visualization 
of q uantum information protocols in discrete phase-space 
ElBQEEEElm; further work might bring 
insights into existing applications, or suggest new ones. 
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